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The interrelation between spin and charge in semiconductors leads to interesting effects, e.g., the 
Rashba-type spin-orbit splitting or the exchange enhancement. These properties are proposed to be 
used in applications such as spin transistors or spin qubits. Probing them on the local scale with 
the ultimate spatial resolution of the scanning tunneling microscope addresses their susceptibility 
to disorder directly. Here we review the results obtained on two-dimensional semiconductor systems 
(2DES). We describe the preparation and characterization of an adequate 2DES which can be probed 
by scanning tunneling microscopy. It is shown how the electron density and the disorder within the 
2DES can be tuned and measured. The observed local density of states of weakly and strongly 
disordered systems is discussed in detail. It is shown that the weakly disordered 2DES exhibits 
quantum Hall effect in magnetic field. The corresponding local density of states across a quantum 
Hall transition is mapped showing the development from localized states to extended states and 
back to localized states in real space. Decoupling the 2DES from screening electrons of the bulk 
of the IH-V semiconductor leads to a measurable exchange enhancement of up to 0.7 meV which 
depends on the local spin polarization of the 2DES. At stronger confinement potential, i.e. larger 
doping, the Rashba spin splitting with a as large as 7 ■ 10~^^eVm is observed as a beating in the 
density of states in magnetic field. The Rashba spin splitting varies with position by about ±50 % 
being largest at potential hills. 
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I. INTRODUCTION 

The proposal of the spin transistor by Datta and Das 
in 199C|i' initiated a new field, which is called semicon- 
ductor spintronics. It has the goal to use the spin degree 



of freedom of the electron not only to store and read 
information, but also to process information.'^ The 
initial idea of the spin transistor was to use the electric 
field of a gate electrode to induce an effective magnetic 
field -Bgff within the reference frame of moving electrons 
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of a 2D electron system (Rashba effect)!^ This effective 
B-field induces precession of electron spins, if they are 
oriented perpendicular to B^g and if the electrons move 
perpendicular to 13 ^g. Using ferromagnetic source and 
drain electrodes, which are polarized both in the same 
direction perpendicular to B^^g, the transistor can be 
tuned on or off, if the precession angle between source 
and drain is 27r or tt, respectively. This solely depends 
on the gate voltage and, thus, voltage pulses can be 
used for tuning. Such a spin transistor might be faster 
and more energy efficient than standard field effect 
transistors. 

Motivated by this idea, various experimental tech- 
niques using electron spin resonance ,1^"^ time-resolved 
optics,'^'^ and multi-terminal spin transporifi^Hi^ have 
been applied to un derstand the spin injection into 
semiconductors E^liMIll dynamics and possibilities 
of manipulation within semiconductor d^^"^^ as well 
as the spin detection.-^'-^- While highly efficient spin 
injectiorfSEHlj large spin relaxation lengthd^ as well 
as an efficient manipulation of spins by the Rashba 
effedP^ltSi] have been demonstrated, the realizati on of a 
reliably operating spin transistor remains elusive .'^21211 
Another branch of application in semiconductor spin- 
tronics is qu antum computation based on spins in 
quantum dots.^^^ Using the spin degree of freedom of 
an electron as a qubit of information has the advantage 
that the spin is only marginally coupled to other degrees 
of freedom mostly via the spin-orbit interaction, the 
hyperfine interaction to the nuclear spins and the dipolar 
coupling to other electron spins. Thus, decoherence and 
relaxation of the qubit is weak, which might allow to 
implement adequ ate er ror correction schemes for infor- 
mation processing .1221311 Using quantum dots as hosts for 
individual electron spins allows to transfer the immense 
knowledge on preparational and electrical control of the 
quantum dots gained in previous yeaxsW^ However, ade- 
quately manipulating the spins on the local scale remains 
a challenge, since typical electron spin resonance schemes 
are never local on the nm-scalcl^S Thus, researchers 
used singlet and triplet degrees of freedom of a pair of 
electron spins to perform the single qubit operations 
electrically by gates, thereby implementing a promising, 
scalable scheme .l^^l Coherence times up to 0.2 ms have 
been reached, which would lead to about 10^ possible 
qubit operations within the coherence time.'^S! The 
fidelity of a two-qubit operation has been determined 
to be 70 % so far and, thus, still has to be improved.'^ 
Although other approaches to solid-state based qubits 
like superconducting transmon^33MI or NV centers 
in diamoncP^I^ are ahead with respect to fidelity or 
coherence times, it is still likely that the easier scalability 
of the spin qubits renders it the more favorable approach. 

Avoiding disorder in electronic or spintronic devices 
is typically advantageous, although exceptions like the 
quantum Hall resistance standard'*'^ relying on a remain- 
ing amount of disorder are known. More importantly, dis- 



order can never be avoided completely and, thus, the re- 
lation of the relevant effects to disorder are important to 
be understood. Since disorder is a local property, getting 
knowledge on the local scale is the most straightforward 
approach. The highest spacial resolution is provided by 
the scanning tunneling microscope (STM), SI which in dif- 
ferent spectroscopic modes is sensitive to the charge and 
spin distribution of individual electronic states.-^*^ For 
certain structures, even the phase of the electron wave 
function has been reconstructed.'^ A quasi-coherent ma- 
nipulation of a pair of wave functions within a quantum 
corral has also been realized by manipulating the local 
disorder via STM.l^ Thus, STM is an excellent tool to 
probe the interaction of relevant spin properties with dis- 
order down to the atomic length scale. 
To this end, we will show that the exchange interac- 
tion between spins as well as the Rashba-type inter- 
action of spins with electric fields can be detected in 
semiconductors on the nm-scale. Therefore, we use a 
two-dimensional electron system (2DES) which is in- 
duced by surface doping of a low band gap HI-V 
semiconductor.'SsHIII This 2DES is located directly below 
the surface (depth: 10 nm) and, thus, can be probed 
by photoelectron spectroscopj^ as well as by STIVP^HH], 
This way, the parabolic dispersion of the 2DES can be 
measured directly. Moreover, the potential disorder can 
be determined with meV and nm accuracy using the so- 
called tip-induced quantum dot.^^The 2DE S exhi bits the 
integer quantum Hall effect in magnetic fielcP^'^ and the 
corresponding states across a quantum Hall transition 
have been detected by STM.I^Sl Since the magnetic field 
allows to localize electrons, a sizable exchange interaction 
between the localized spins appears, which can be deter- 
mined by careful inspection of the spin splitting within 
local Landau level fans.l^ The spin splitting is roughly 
0.7 meV larger, if the local filling factor is odd than if 
it is even. Moreover, by preparing the 2DES such that 
it is confined within a large electric field (2 • 10^ V/m), 
we can detect the Rashba-type spin splitting. This spin 
splitting is visible as a beating pattern of the density of 
states (DOS) in magnetic field, which without Rashba- 
interaction exhibits regularly spaced Landau levels.l^ 
The Rashba-type spin splitting can also be probed di- 
rectly as a locally fluctuating energy distance between 
two spin levels appearing as peaks in the local density 
of states.lsni Thus, we demonstrate that both, exchange 
interaction as well as Rashba-type interaction of electron 
spins can be probed down to the length scale of individ- 
ual electronic states. 

This review is organized as follows. Chapter [TT] de- 
scribes the characterization of the 2DES by photoemis- 
sion. Chapter |HI| is a short general introduction into 
scanning tunneling spectroscopy (STS). Then, techniques 
to determine the potential disorder of the 2DES with high 



resolution are introduced (chapter IV ) . The resulting lo- 



cal density of states reacting onto different strength of 
the disorder differently is shown in chapter |Vj After de- 
scribing the observed local density of states within the 
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quantum Hall regime (chapter VI), we finally show how with 



the local exchange interaction can be directly extracted 
from the STS data of the 2DES (chapter [Vll| as well 
as how the Rashba effect is probed quantitatively on the 
local scale (chapter VIII). 



II. PROBING THE 2DES BY PHOTOEMISSION 

The general problem of studying 2DESs by surface sci- 
ence techniques is, that they are typically located about 
100 nm below the surface. This depth is chosen inten- 
tionally in order to avoid scattering at surface defects 
and, indeed, high mobilities up to 10'* m^/Vs of the 
corresponding 2DESs have been achieved at low tem- 
perature implying an electron mean free path close to 
mm.'^SHU] However, these 2DESs can not be tackled by 
surface science techniques as STM or angular resolved 
photoelectron spectroscopy (ARPES). If one estimates, 
e.g., the tunneling current from a tip above the surface 
into a GaAs 2DES covered by 100 nm of AlAs with the 
triangular barrier approximation,'^ one ends up with a 
tunneling current of the order of 10"^*^ A, which corre- 
sponds to about 10"^'^ electrons within the anticipated 
age of the universe, well below the detection limit of avail- 
able current amplifiers. This impossibility is illustrated 
in Fig. [ija). One possibility to overcome the problem 
is to tunnel into the cleaved edge of the samp le, wh ich 
gives access to the cross section of the 2DES.'^^'^ In- 
deed, subbands and confined states of the 2DES have 
been mapped by STS of cleaved heterostructures.I^S] Hov- 
ever, the distribution of electron wave functions within 
the 2DES plane, guiding, e.g., the localization phenom- 
ena of the 2DES, is not visible by this method. Another 
possibility is to induce a 2DE S direct ly at the surface by 
surface doping (see Fig. [lJb)).'^^'^^^This is possible for 
the low bandgap materials as InAs and InSb, which ex- 
hibit a relatively large electron affinity.^^ Basically, most 
adsorbates deposited on that surfaces exhibit a confined, 
occupied st ate ab ove the conduction band edge of the 
substrate.' ^^^ * ^^! '^ Thus, the electron is moved from 
the adsorbate into the semiconductor leaving a positive 
charge behind, which induces confinement of electronic 
states in the conduction band. The resulting band bend- 
ing in the near surface region together with the wave 
function of the confined state with lowest energy (^1(2:) 
of the first subband) is sketched in Fig. ^c). The band 
bending as well as the wave function can be calculated 
using so-called Poisson-Schrodinger solvers®^, which take 
the screening of the charged surface adsorbates by bulk 
dopants as well as by confined electrons self-consistently 
into account. A reasonable approximation is given by 
assuming a triangular potential well in z direction with 
infinite surface barrier leading to'^ 



6 = 



48 • 7rm*me • • N ^^^^ 
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N = A^Dopant + T^A^2DES- 



(2) 
(3) 



*(z) = (^)0-5 



z ■ e 



-bz/2 



(1) 



Here, m* is the effective mass of the conduction band in 
^-direction, e the dielectric constant of the semiconduc- 
tor, iVDopant the 2D density of charged dopants within 
the depletion or accumulation layer and A^2DES the 2D 
density of confined electrons (eq = 8.8 • 10~^^ As/Vm, 
n = 1.05 • 10-34 Js, me = 9.1 • 10-31 kg). Typically, the 
center of mass of ^'(z), zq = 3/6, is about 10 nm below 
the surface and, more importantly, the 2D electrons spill 
out of the surface due to the finite barrier height at the 
surface. Thus, they can be detected by STM like surface 
states, e.g. on metals.'^ 

A first characterization of the induced band bending and 
the resulting 2DES subbands can be performed using 
ARPES. Fig. [ijd) shows a number of spectra obtained 
for different coverages of Nb on an n-type InAs(llO) sur- 
face. Thereby, InAs(llO) is used since other surfaces as 
(111) or (001) typically exhibit a band bending already 
after surface preparation without adsorbates,'S2i presum- 
ably because of the relatively large defect density. In 
contrast, the InAs(llO) surface exhibits only step edges 
as natural defects with typical interstep distance of 5 
/im.^3 Thus, it shows flat band conditions after cleavage 
in ultrahigh vacuum. This is shown, e.g. in Fig. [ijg), 
where the valence band of p-type InSb(llO) as measured 
by ARPES exhibits heavy and light hole bands as well 
as a spin-orbit split band (maximum at —0.8 eV) and 
a Fermi level E-p that cuts the valence band, which evi- 
dences the absence of band bending. Consequently, ad- 
sorbates can be used to tune the band bending. The 
subsequently deposited adsorbate density is determined 
by a quartz balance and is given with respect to the unit 
cell of InAs(llO), i.e. 1 % coverage corresponds to an 
adsorbate density of 3 • 10^^ m"^. The three spectra in 
Fig. [l|d) all exhibit two peaks, which can be identi- 
fied as bulk valence band states of InAs.f^ They shift to 
lower binding energy with increasing coverage indicating 
the adsorbate induced band shift in the near surface re- 
gion. Fig. [ije) shows the binding energies of a number 
of ARPES peaks observed at different photon energies 
as a function of Nb coverage. Generally, all peaks show 
similar energy shifts. However, the ones with the largest 
escape depths of the photoelectrons show an about 20 % 
smaller peak shift in line with the idea that they probe 
the most shifted bands deeper into the bulk of the sys- 
tem. The fact that the measured work function shift is 
larger than the observed peak shifts, in particular, at 
large coverage implies an additional change of a surface 
dipole, which has been attributed to the lifting of the re- 
laxation at the InAs(llO) surface.'^ In the case of Nb on 
n-type InAs(llO), the largest band bending is 320 meV. 
It is most likely limited by the donor level of the adsor- 
bate being i?ads = 320 meV above the conduction band 
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FIG. 1. (a) Sketch of buried 2DES as typically used in transport measurements with mobility /i and calculated tunneling 
current from the tip I indicated; (b) same as (a), but for an adsorbate induced 2DES; (c) calculated band structure of InSb 
after depositing Cs atoms on top; Fermi level Ep, charged Cs on the surface, subband energies (horizontal lines) and wave 
function 'i'{z) of the first subband are shown; the valence band is omitted for clarity; (d) ARPES data of n-InAs(llO) (energy 
distribution curves) {Nd = 1.1 • 10^^ m"^) with different Nb coverages on the surface as indicated; nearly vertical lines mark 
the shifting peaks, photon energy hv = 13 eV, emission angle of photoelectrons = 0°; (e) energy shifts of different valence 
band peaks of n-InAs(llO) as a function of Nb coverage, lower right inset shows photon energy hv, resulting kinetic energy i5kin 
of the photoelectrons and estimated escape depths A of the photoelectrons, the surface dipole determined by the differences 
between peak shift and workfunction shift is marked, notice that the peaks shift to lower binding energies with increasing 
coverage; (f) highest position of the Fermi level above the conduction band minimum at the surface for different adsorbates 
on InAs(llO) as a function of ionization energy of the free adsorbate atomf^ (g) ARPES data of p-InSb(llO) {Na = 1.5 • lO^"* 
m~^) directly after cleavage; (h) same data as (g) after adsorbing 2 % Cs, T = 80 K, hv = 21 eV. 



minimum of InAs. If the band bending is large enough 
such that the donor level is shifted below the Fermi level 
of the substrate, no additional charging of the adsorbates 
occurs. Thus, the process is self limiting. In turn, the 
adsorbate level can be determined via the largest band 
shift at the surface. 

It has been found that this adsorbate level scales inversely 
with the ionization energy of the corresponding neutral 
atom as shown in Fig. n|f).'5^ Obviously surface band 
shifts up to 600 meV are possible using alkali atoms as 
adsorbates. This is much larger than the band gap of 
InAs being 430 meV at T = K or 360 meV at T = 300 
K.ISS Such a large band bending can even induce occu- 
pied confined states within a p-doped InAs crystal, i.e. 
a conducting 2DES within an inversion layer. Fig. ^g) 
and (h) show a direct proof of such an inversion layer 
using ARPES mapping of the E{k) dispersion of p-type 
InSb(llO) prior and after deposition of 2 % of Cs. The 
valence band of InSb is visible in Fig. [ijg) as described 
above. Adsorbing 2 % of Cs moves the valence band 
down by about 450 meV and a single parabolic band 
curved upwards in the E{k) dispersion appears at Ep 
(Fig. This band is considerably weaker in inten- 



sity than the valence bands, such that the valence band 
appears as a completely dark area, if the contrast is cho- 
sen such that the 2DES is visible. This is due to the 
fact that ARPES probes only a few percent of the 2DES 
charge density, since the photoelectron escape length is a 
factor of 10-20 smaller than the extension of the 2DES. 
The fact that the energy distance between the origin of 
the upwards parabola of the 2DES to the valence band 
maximum is 350 meV, i.e. much larger than the known 
band gap of InSb (230 meVJP^ directly reflects the con- 
finement energy of the first subband of the 2DES being 
120 meV. Higher subbands of this 2DES are not occupied 
and, thus, not visible in ARPES. 

Importantly, a high level of control of the band struc- 
ture properties of the 2DES is obtained by ARPES. The 
2DES properties can be designed by varying the dopant 
type (n- or p-type) and dopant density of the substrate 
as well as the adsorbate type and coverage. 
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III. SCANNING TUNNELING MICROSCOPY 

Before describing the STS data of the 2DES, we will 
shortly summarize the most important features of STS. 
In STS, a sharp metallic tip, mostly ending in a single 
atom, is positioned 3 — 8 A above a conducting surface. 
The surface is usually prepared in ultra high vacuum 
(UHV) in order to be atomically smooth. A voltage V is 
applied between the tip and the conducting surface and 
the resulting tunneling current / is measured. / depends 
exponentially on the distance between surface and tip Az 
according to I{z) oc e^"^^ . A good approximation for a 
is a = A/(4mc • ($s + <&t - e|y|)/7i) with Planck's con- 
stant h, electron mass rric, electron charge e and work 
functions of tip and sample $f and <i>s, respectively. A 
good estimate is a ~ 2.1/A. 

The tip is positioned with respect to the sample using 
piezoelectric elements. All three directions x, y, and z 
can be changed with sub-pm precision.!^ For STM, the 
tip is scanned in x and y direction and the tunneling cur- 
rent is kept constant by a feedback mechanism adjusting 
z. The resulting z{x,y) is plotted and called constant- 
current image. It represents, to first order, a contour of 
constant integrated local density of states of the substrate 
LDOS{x,y, z, E), where the integration has to be taken 
between the Fermi levels of sample and tip to be adjusted 
by y.l^The central position of the very last atom of the 
tip is given by {x, y, z) and E is the energy. Such images 
are often called topography of the sample indicating that 
corrugations of the atomic positions dominate the con- 
tour. 

Differentiating / with respect to V (at low V with respect 
to ^s{t)) gives direct access to the LDOS according to: 



dI/dV{x, y, z, V) oc LDOS{x, y, z, E) = 
Y,\^'^{x,y,z)\'-5{E~E). 



(4) 



Thereby, ^' ^ are the single-particle wave functions of the 
substrate at energy E and E = e ■ V . Of course, this 
requires that the system is adequately described by in- 
dependent single-particle wave functions. Moreover, an 
s-type symmetry of the orbitals of the last atom of the tip 
is necessary, in principle. Chen has shown that the model 
remains largely correct even if higher orbital momenta of 
the tip wave functions are contributing.!^ In real exper- 
iments, the (5-function has to be replaced by an energy 
resolution function with approximate full width at half- 
maximum of 5E « ^(3.3 • keTf + (1.8 • eKnod)^- T is 
the temperature {k^: Boltzmann's constant) and Vmod is 
a modulation amplitude used to detect dl/dV by lock- in 
technique. The resulting dl /dV{x,y) recorded at con- 
stant V and constant / is often called LDOS-im&ge im- 
plying that the image is proportional to the LDOS distri- 
bution in the surface layer. Strictly speaking, this is not 
correct since constant / is not identical to constant height 
above the surface z{x,y) ^ const.. A local enhancement 
of density of states, e.g. around a dopant within a semi- 
conductor, can locally lead to a larger tip-surface distance 



and, thus, to a reduced differental conductivity. This can 
be cured by measuring the global distance dependence of 
the tunneling current I{z) and dividing dI/dV{x,y) by 
I{z{x,y)) resulting in: 



LDOS{E,x,y) 



dI/dV{x,y,V = E/e) 
I{z{x,y)) 



(5) 



However, that requires that the decay of / with respect to 
z is spacially constant^^, which is a reasonable assump- 
tion for the LDOS of a 2DES, but could fail on the atomic 
scaleP. 

Another modification, which applies favorably for 
dl /dV{V) curves at higher V with respect to ^s(t)^ is 



LDOS{E) cx 



dI/dV{V ^ E/e) 

WW 



(6) 



This compensates for the effect that the transmission 
coefficient of electrons through the tunnel barrier also 
changes with applied voltage, which is relevant if |V^| > 
0.5 V.'SSI Another problem with dI/dV{V) curves is that 
they are sensitive to both, the LDOS{E) of the sample 
as well as the LDOS{E) of the tip. Thus, it is required 
to have a featureless LDOS{E) of the tip. Typically, 
one repeats measurements with different tips on the same 
sample and attributes repeatedly observed features to the 
sample or one changes the properties of the sample, e.g. 
by applying a -B-field and measure them with the metal- 
lic tip being insensitive to the _B-field.'^ 
Thus, STM can measure atomic structure with sub-pm 
resolution and electronic structure (LDOS) with sub- 
meV resolution.lSSl The energy resolution makes STM 
complementary to the transmission electron microscope 
(TEM) which reveals atomically well defined structural 
information partly with chemical specificity,!^ but not 
the LDOS down to the meV scale. 



IV. DISORDER POTENTIAL 

A. Experimental determination 

One drawback of STS on semiconductors is the low 
charge density of semiconductors and the resulting rela- 
tively large screening length Ag of several nm, which can 
be calculated according to: 



ne ■ Ob 



Am* 



(7) 



with the dielectric constant e = 14.6 for InAs and 
e ~ 16.8 for InSb, the Bohr radius as = 0.053 nm, 
and N being the 3D density of mobile charge carriers. 
Thus, electric fields from the tip can penetrate into the 
semiconductor leading to the so-called tip-induced band 
bending. ^■^ 

This band bending depends on the work function differ- 
ence between tip and sample A<I> and the applied voltage. 
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FIG. 2. (a) Spatially averaged dl/dV curve of n-InAs(llO) (iVn = 1.1 ■ 10^^ m"^) after depositing 2.7 % of Fe, Fermi level 
-Ef and the subband energies E\, E2 as deduced from ARPES data are marked by vertical lines; (b) spatially averaged dl/dV 
curves of the same n-InAs(llO) as in (a) probed with a different tip prior (lower curve) and after (upper curve) depositing 4.5 % 
Fe, features belonging to the 2DES, the 3DES and the tip induced quantum dot (QD) are marked; (c) grey scale plot of dl/dV 
as a function of voltage V and position x, the subband energies Ei, E2 and the feature belonging to the QD are marked, the 
latter is used to map the electrostatic potential of the 2DES; (d), (e) measured electrostatic potential landscape of the 2DES 
induced by 0.8 % and 4.5 % Fe, respectively, greysacle covers a range of 20 meV from black to white, note the larger disorder 
potential after depositing 0.8 % Fe; (f) STM image of the area marked in (e) (1/ = 0.1 V, 7 = 50 pA), the black dots are the 
positions of Fe atoms; T = 6 



Typically, the work function of the tip varies by about 
300-500 meV depending on details of the microtip.'^ 
Thus, it is important to prepare the tip by trial and error 
such that A(f> ~ meV. 

On the other hand, the tip-induced band bending can 
be used to determine the disorder potential of the 2DES. 
Therefore, one has to realize that the band bending, if 
downwards towards the surface, leads to a conduction 
band region, which is confined in all three directions, 
i.e. only the bands below the tip are pushed downwards 
by the electric field of the tip leading to a quantum dot 
below the tip. This tip-induce d qua ntum dot (TIQD) 
can be moved across the suriaceW^ In z-direction per- 
pendicular to the surface, the confinement region of the 
quantum dot states is approximately the same as the 



confinement region of the 2DES, since screened by the 
same dopant density. The corresponding probability am- 
plitude is given by eq. [ij however with slightly differ- 
ent A^2DES- In (a;, y)-direction parallel to the surface, 
the confinement is less, since the tip has a certain lat- 
eral extension, which adds to twice the screening length. 
A good approximation is a Gaussian with a full width 
at half maximum (FWHM) of 50-100 nm. This width 
can be determined by the energy distance of the con- 
fined states which appear as peaks in dl/dV curves. The 
lowest-energy state of such a quantum dot has a Gaussian 
lateral shape according to 



|vI/i(x,y)|2oce (8) 
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with a FWHM 2 • cti about 4-5 times smaller than the 
FWHM of the QD, i.e. about 10 nm.l^^ 
If the TIQD is absent, dl/dV curves measure the unper- 
turbed LDOS of the sample. An example of the spatially 
averaged dl/dV curve of an adsorbate induced 2DES 
without TIQD is shown in Fig. exhibiting two steps 
according to the two low-energy subbands of the 2DES. 
Above the steps, the spatially averaged LDOS is flat as 
expected from the density of states of each subband of a 
2DES 

„ , , TTlp • TO* , , 

Dn{E) = ^2 ■ (9) 

The second subband appears less intense, since it is 
buried deeper into the substrate. Importantly, the 2DES 
leads to dl/dV intensity in the region below Ep although 
the InAs-sample is n-doped, i.e. the 2DES is within the 
region where the bulk band gap of the semiconductor is 
expected. 

In case of a TIQD, already 3DES before the surface dop- 
ing by adsorption exhibits peaks below Ep, i.e. in the 
band gap region of the 2DES. From the distance of the 
peaks of about AE = 40 meV, one can calculate the 
width of the quantum dot and the FWHM of the lowest 
state using the harmonic oscillator approximation 

AE = huj ^ (7i = ^= ^ ^ =8.5nm, (10) 
V'Tie • m* ■ AE 

i.e. the lowest state has a FWHM of 17 nm. 
Moving this state across the surface exposes it to the 
potential disorder of the sample. Thus, the confinement 
potential of the QD is modified by the 2D disorder poten- 
tial of the 2DES V{x,y) = J 'i/i{z)V{x,y, z)dz , where 
V{x,y,z) is the bare 3D electrostatic potential of the 
sample. Thus, within first order perturbation theory, the 
energy of the lowest state is 

Ei{x,y) =< 'i>i{x-x,y-y)\V{x,y)\-^i{x-x,y-y) 

(11) 

Consequently, an energy map of the lowest state is a map 
of the electrostatic potential of the 2DES, which is coars- 
ened by 2(7. Notice that the 3D potential is already coars- 
ened by ~ 3/& ~ 10 nm due to the extension of the 2DES 
in z direction. Thus assuming an isotropic screening of 
point like charges as the origin of the disorder, the lowest 
state maps the potential down to about the lowest rele- 
vant length scale. 

Fig. ^c) shows a dI/dV{V, x) map of an InAs-2DES as 
a greyscale plot. The subband energies i?„ determined 
from ARPES of the same system are marked. At volt- 
ages above these energies, the 2DES shows wavy patterns 
along the lateral position x which are related to the scat- 
tering of the electron waves (see chapter |v]) . Since more 
states contribute to the scattering above E2, the apparent 
wave length gets smaller. Below Ei a strong peak mean- 
dering in energy as a function of x is visible. This peak 
marked QD is the lowest state of the quantum dot (see 
also upper curve in Fig. ^b)). In the central x region at 



70-130 nm even the second QD state is visible showing a 
similar meandering as the first one. If one plots the peak 
voltage of the low energy state as a function of position, 
one gets the E{x, y) maps shown in Fig. [2jd) and (e) for 
two different 2DES, respectively. These images are called 
map of the potential disorder. They are both measured 
on n-InAs with a donor density of iV^ = 1.1 • 10^^ m~'^, 
but with different densities of Fe adsorbates as marked. 
The FWHM of the histograms of all potential values is 
20 meV in both cases, but surprisingly the potential dis- 
order varies on a shorter length scale at lower adsorbate 
coverage. 

Fig. [2]je) features about 4 deep troughs and about 7 shal- 
lower troughs. The number of donors within the 2DES 
region of the imaged area A = 200 x 200 nm^ amounts 
to A'd2D = A ■ N-£, ■ 3/b = 4 implying that the deep 
troughs mark the donors of the substrate located within 
the 2DES. The shallower troughs are most likely sub- 
strate donors lying below the 2DES, which are not com- 
pletely screened. This is compatible with As = 20 nm for 
the bulk of InAs being a factor of two larger than 3/6, 
although As is slightly shortened by the presence of the 
2DES. 

Thus, one concludes that the surface dopants are not rel- 
evant for the potential disorder. This gets even more 
obvious, if one compares the distribution of Fe atoms, 
which are visible as black dots in Fig. ^i). The poten- 
tial fluctuations appear on a length scale covering several 
100 adsorbate atoms. More precisely the 2D density of 
dopants within the 2DES is a factor of 1000 lower than 
the density of Fe atoms. The potential troughs in Fig. 
^d) are about 80/(400x400) nm"^, i.e., at least, afactor 
of two denser than in Fig. ^e). Thus, they are obviously 
influenced by the adsorbate layer. The likely reason for 
this different behavior will be discussed below. 
Importantly, we can probe the disorder potential experi- 
mentally on the meV and nm scale.'^ Unfortunately, that 
requires the presence of a disturbing TIQD. Thus, the 
pristine LDOS of a 2DES with exactly known potential 
disorder cannot be measured. 

Nevertheless, we mapped the disturbed LDOS at differ- 
ent energies and compared it with the calculated LDOS 
resulting from a diagonalization of the matrix 

<^^i^'y)\^—^ + V{x,y)\'^j{x,x) >, (12) 

where Vl'i(j) are plane waves, and we assumed periodic 
boundary conditions. In order to get the LDOS, the re- 
sulting wave functions are averaged using a Gaussian en- 
ergy distribution with FWHM= SE . The length scale of 
LDOS fluctuations, the histogram of LDOS values (see. 
Fig. [5jc)) and the Fourier transformation of the mea- 
sured LDOS are well reproduced, but only few of the 
exact shapes of the LDOS(x, y) are identical leading to a 
cross correlation of 15 %^ Possible reasons are the re- 
maining influence of the TIQD, the coarsening of some 
of the details of the potential disorder and a localization 
length of the states which might be larger than the im- 
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age size. Indeed, taking the measured mean free path 
of 300 nm (see below) , an influence of the last argument 
appears likely. 



B. Role of surface dopants 

Next, we will discuss the role of the adsorbates on the 
disorder. Therefore, it is essential to realize that only 
part of the adsorbates are charged. The requirement for 
charging is that the electric potential produced by all 
other adsorbates at the position of a given one is lower 
than i^ads — ^^BCBM- Neglecting the screening in the ad- 
sorbate layer, one can easily simulate the charging prob- 
ability of an adsorbate using £^ads — -Ebcbm — 320 meV 
and Xs = 10 nm by depositing adsorbates subsequently 
onto random positions a;„ on a 2D plane and surrounding 
them by 
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n-l 

E 



E 



BCBM 



(13) 



(14) 



and Vn{x) = eV otherwise. One has to take into ac- 
count that e at the surface is reduced with respect to the 
bulk.l^ The result for a random choice of a;„ is shown in 
Fig. |3][a) as black triangles. For comparison the ion- 
ization probability for Fe-adsorbates on n-InAs(llO)^ 
is shown as open circles revealing reasonable agreement. 
It has been calculated from the measured band bending 
by solving the Poisson-Schrodinger equation and dividing 
the sum of the resulting 2DES density A''2des , which has 
also been determined experimentally,'^ and the density 
of neutralized donors iVDopant by the adsorbate density. 
Importantly, the ionization probability is 40-50 % at 0.8 
% coverage but only 8 % at 4.5 % coverage. Thus, the 
2DES density changes only slightly by increasing the cov- 
erage to 4.5 %, while most of the additional Fe atoms are 
uncharged. The remaining electrons in the Fe layer can 
screen the charged Fe-atoms effectively. The charge den- 
sity in the Fe layer of about 1.5 • 10^^ m~^ can probably 
even rearrange dynamically, since it is energetically cou- 
pled to the 2DES. Given the very effective and density 
independent screening in 2D systems above a threshold 
density,'^ it is reasonable to assume that the charge den- 
sity in the Fe layer is effectively screened on the length 
scale of the average Fe-Fe distance being 2.5 nm at 4.5 
% coverage. Thus, the 2DES which is coarsening on the 
length scale of 10 nm due to its extension perpendicular 
to the surface basically ignores the potential disorder of 
the Fe layer provided that the adsorbate density is sig- 
nificantly larger than the density of transferred charge. 
This simple picture breaks down, if the adsorbates are 
mobile on the surface and, thus, form clusters. Basically, 



each cluster is only charged once due to Coulomb block- 
ade. Thus, the ionization probability per adsorbate is 
reduced, if clusters are formed and the maximum band 
bending is achieved only at higher coverage. Moreover 
the ionization level could be changed by hybridization. 
The squares in Fig. |3][a) show the ionization probability 
for Nb which forms cluster consisting of 2-4 atoms de- 
pending on coverage after room temperature deposition. 
The complex curve deviating strongly from the monomer 
model implies that both, single charging of multimers as 
well as hybridization in multimers is essential.'^ 
Importantly, already initial clustering of adsorbates leads 
to additional charge inhomogeneities. This is due the 
fact that cluster formation implies a reduced adsorbate 
density around the cluster and, thus, leads to a local 
potential hill due to reduced charging. This has been 
experimentally evidenced by Hashimoto et al.W^ Thus, 
it is important to avoid clustering and to deposit sig- 
nificantly more adsorbates than required for the charge 
transfer into the 2DES, if one wants a high mobility 
2DES. Measured values of the Hall mobility ii after 
Ag deposition on p-type InSb(llO) {Na = 1 - 2 • lO^^ 
m~^) are up to ^ = 5 m^/Vs at a 2D carrier density 
iV2D = 1-5 • 10^^ m~^^ This implies a mean free path 
Amfp of Amfp = ^/e ■ •v/27r7V2D ■ M — 320 nm. Hall mea- 
surements at 40 % Fe coverage or 1-20 % Ag coverage 
on top of higher p-doped InAs(llO) {Na = 1.7 • lO^^ 
m^"^) revealed a lower mobility of 0.7 m^/Vs and 0.3 
m^/Vs, respectively, again pointing to the leading im- 
portance of bulk doping for the disorder potential. ^^-^ 
In accordance with the above analysis, it was found that 
the mobility decreased considerably below a Ag coverage 
of 1.5 %. Interestingly, spin-glass behavior in the Fe/p- 
InAs(llO) system has been deduced from hysteresis and a 
logarithmic time dependence of relaxation dynamics be- 
ing present up to about 1000 s after a change of B field. 
This effect has been correlated with the complex inter- 
play of antiferromagne tic an d ferromagnetic coupling in 
Fe-chains on InAs(llO) .123211 

Fig. [31(c) shows the ionization probability for Co on p- 
InAs(110), which forms larger islands at room tempera- 
ture as shown in the inset. The ionization probability is 
plotted as a function of the island density and evidently 
gives a very similar curve as the one found for the Fe 
adsorbate density in Fig. [3](a). Thus, indeed, also larger 
islands are only charged once. The reason is the Coulomb 
blockade which can be directly seen in dI/dV{V) curves 
recorded on a particular Co island as a gap surrounded 
by two peaks (see Fig. [3|^c)).''°. A histogram of the size 
of the Coulomb gaps has been performed for an island 
density of 3.5 • 10^^ m~^, where each island consists on 
average of 15 atoms leading to a an average gap size 
of 465 meV with a standard deviation of 150 meV. An 
estimate using a metallic sphere with the radius of the 
Co-clusters i? in a distance d from a metallic plate re- 
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FIG. 3. (a) Open symbols: ionization probability per adsorbate atom for Fe and Nb deposited on n-InAs(llO) (A'd ~ 1-1 • 10^^ 
m~^) as deduced from the coverage dependence of the band shift measured by ARPES (see Fig. [l^), full symbols: ionization 
probability calculated from a random distribution of monomers with an ionization level 300 meV above the conduction band 
minimum of InAs, inset: STM image of Fe adsorbates on InAs(llO) (/ = 200 pA, V = mV); (b) Ionization probability as 
a function of Co island density on p-InAs(llO) (A'^a = 4.6 • 10^"^ ni"'^), inset: STM image of Co-islands on InAs(llO) prepared 
at T = 120° C (7 = 50 pA, V = —0.86 V); (c) dl/dV spectrum obtained on a single Co island exhibiting a Coulomb gap 
(Vstab = 0.7 V, /stab = 0.7 uA, Kiod = 5 mV); (d) simulated potential landscape of p-InAs(llO) {Na = 4.6 • 10^^ m"^) covered 
with 15 % Co atoms per InAs(llO) unit cell, respectively, a Co island density of A^isiand = 3.5 • lO'^^ m"^, it is assumed that 
only the 50 % of the islands with the lower Coulomb gap are charged; (e), (f) dl/dV curves obtained from the areas marked 
in (d); (\4tab = 0.7 V, /stab = 0.7 nA, ymod = 5 mV); subband energies Ei, E2, E3 of the 2DES and the peaks of the strongly 
localized states (Is, Ip) are marked; (g), (h) dl/dV im ages t aken at the peak energies of the localized states marked in (f) 
(Fstab = 0.43 V, /stab = 0.7 nA, V^^i = 5 mV); T = 6 K.™ 



InAs of Ai?chargc = fEpf ' g:^^ ^ The calculated values 
are about a factor of two larger than the value of the 
measured Coulom b g ap. This might be caused by the 
background chargeiSSI as observed, e.g. for Co clusters on 
MgO(0001),'i^ or by an extension of the wave functions 
of the metal cluster into the InAs. A crude estimate of 
the 2D disorder potential can be made by assuming that 
all the clusters with the lower Coulomb gap up to the ion- 
ization density are charged and, are, thus, surrounded by 
a screened Coulomb potential. Adding charged acceptors 
randomly and averaging as usual over the z-extension of 
the 2DES leads to the 2D disorder potential shown in 
Fig. ^d). It fluctuates by about 150 meV on a length 
scale of 20 nm. Thus, the potential fluctuation is an or- 
der of magnitude stronger than in Fig. ^e) . This allows 
to probe the 2DES at very different disorder strengths. 



V. LOCAL DENSITY OF STATES AT 
DIFFERENT STRENGTH OF DISORDER 

The two-dimensional electronic system (2DES) at the 
surface induced by adsorbates shows a nearly parabolic 
dispersion: 

E{k)-E^^,^^ (15) 

with E{k) being the energy of the state. En the subband 
energy corresponding to the confinement vertically to the 
surface, k the wave vector of the electron parallel to the 
surface, uie = 9.1 • 10~^^ kg the bare electron mass, and 
m* the effective mass of the conduction band. Thus, 
the system is a good paradigm for free electrons, which 
also show parabolic dispersion and exhibit an extremely 
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0.8 % Fe/n-lnAs(11Q) 




1.5 % Cs/p-lnSb(110) 

FIG. 4. (a)-(d) dl/dV images at the voltages indicated of n-InAs(llO) (TVd = 1.1 • 10^^ m"^) covered with 0.8 % Fe (V^tab = 0.1 
V, /stab = 0.3 nA, Vniod = 1-8 HiV) , insets: Fourier transformation of the real space data with preferential \k\ vector marked in 
(d); (e) potential landscape of the same area determined by the lowest tip induced quantum dot state as shown in Fig. [2|c)-(e); 
(f) preferential |A;| values as deduced from the Fourier transformations in comparison with unperturbed InAs band structure as 
calculated within k ■ p-approximation; (g) dl/dV image at V = V of p-InSb(llO) {Na — 1 — 2 ■ 10^* m""') covered with 1.5 
% Cs (Vstab = 0.3 V, /stab = 0.2 uA, Vmod = 1-0 mV), lower inset: Fourier transformation of the real space data, upper inset: 
STM image of the same surface (/ — 30 pA, V = 0.3 V) shown to scale and proving that the dark dots in the dl /dV image 
are caused by the Cs atoms; (h) ARPES data of p-InSb(llO) {Na = 1 - 2 • 10^* m"^) covered with 1.5 % Cs, hv = 21 eV; a 
parabolic approximation of the band of InSb with effective mass marked (dash ed lin e) and the preferential |fc|-values from the 
Fourier transformation of STS data are included; T = 6 K (a-g); T = 80 K (h).1^2168l 



interesting phase diagram as a func tion of electron den- 
sity, temperature and magnetic field .1^221 At the energy of 
the bulk conduction band minimum /^bcbMj the effective 
mass of InAs is m* = 0.023 and of InSb is m* = 0.0135.ES1 
In fact, the bands are not perfectly parabolic, but exhibit 
a reduced curvature at higher energy, which can beap- 
proximated by an energy dependent effective mas^^^ 



mn 



1 + 2 ■ {E — /i'bcbm) 



Gap 



(16) 



with the band gap being E'cap ~ 0.43 cV for InAs and 
Eg^p = 0.23 eV for InSb at T = K.'^ 
For a 2DES, this leads toP^ 



TO* (E) — rriQ 



1 + 2.(1/3(K„-^: 



BCBM ) 



E, 



Gap 



, (17) 



i.e. the subband energy perpendicular to the surface E^ 
counts differently than the in-plane energy E^ | . 
The spin-orbit coupling within the valence band couples 
to the conductio n b and, which has preferentially a 5s- 
orbital character .li^SI This can be described, e.g., by k ■ 
p-theory. As a consequence, the (/-factor of the 2DES 



is strongly enhanced, negative, i.e. spin moments align 
preferentially antiparallel to B, and energy dependent. 
A good approximation is;^''^ 



90 



m*{E) ' 



(18) 



with 5o — —15 for InAs and — —51 for InSb. More- 
over, the Rashba-type spin splitting within a potential 
gradient given, e.g. by the confinement of the 2DES in 
z direction, is relatively strong as will be discussed in 
chapter [Villi 

One should keep in mind that the given approximations 
are only valid at energies E — Z^bcbm < Eq^^-^. Impor- 
tantly, the LDOS corresponding to this 2DES can be 
measured by STS, since only a small amount of the sur- 
face is covered by adsorbates. In between these adsor- 
bates the 2DES of the substrate can be probed in dl/dV 
images measured within the bulk band gap of the semi- 
conductor. The result at two different strength of the 
disorder potential are discussed next. 
At strong disorder, some of the electronic wave functions 
of the 2DES get strongly localized, i.e. they are com- 
pletely confined within a single valley of the potential 
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disorder. Approximating the potential trough marked 
by a red dot in Fig. [sjjd) as parabolic with an extension 
of r = 15 nm at 75 meV, i.e. 



(19) 



in good agreement. Thus, obviously the wave pattern 
belongs to the adsorbate induced 2DES. Moreover, it is 
clearly visible that the adsorbates only marginally dis- 
turb the image quality of the 2DES appearing either as 
white spikes in Fig. |4][a)-(d) or as black dots in Fig. gg). 



one gets single particle energies En ~ (n + 0.5) • 40 
meV, which is reasonably smaller than the height of 
the potential trough. Thus, confinement of about the 
first two single particle levels is expected and, indeed, 
observed.'sni Fi gure ^i) shows a dl/dV curve exhibiting 
two sharp peaks with an energy distance of 20 meV. A 
dl/dV map at the corresponding energies, which are 
plotted here for another trough, is shown in Fig. [sf^g) 
and (h) exhibiting an elliptical LDOS distribution for 
the lower peak and a lobe-like distribution as known 
for p-states at the upper peak. Thus, the strongly 
localized states are similar to confined states within 
parabolic quantum dots.^^ The quantitative discrepancy 
of the confinement energy is probably caused by the 
nonparabolicity of the potential trough. Notice that the 
black spots within the LDOS wave functions are caused 
by the presence of the Co islands, which are Coulomb 
blocked. 

Next, we discuss the LDOS data obtained at lower disor- 
der. Fig. |4ja)-(d) shows the LDOS images recorded for 
the potential shown in Fig. [4|e) exhibiting fluctuations 
of about 20 meV only. The total intensity in each image 
corresponds to 40 complete electronic states, but since 
the scattering length of individual states is about the 
image size, more states contribute to the LDOS with 
part of its intensity distribution. 

The LDOS images exhibit corrugations decreasing in 
length scale with increasing voltage. The corrugation 
patterns are rather complicated and do not exhibit 
the simple circular ring structures found in the InAs 
3DES^® or foun d for p arabolic surface states on the 
surface of metals.l^ "'' * ^"^* Fourier transforms (FT's) of the 
LDOS (insets) reveal the distribution of contributing 
fc-values. At low voltages a circle is visible, which 
at higher voltages is confined by a ring. A plot of 
the fc-values corresponding to the rings and, thus, 
dominating the spectrum is shown as a function of 
voltage in Fig. |4]jf). At low voltages where the ring is 
not apparent the outer diameter of the circle is taken. 
For comparison the i<^(fc)-dispersion of InAs according 
to k ■ p theory (see eq. [15] and 17) is drawn which 
exhibits good agreement with the data. However, we 
do not observe a ring only, which is due to the mixing 
of states by the disorder potential according to eq. [T2| 
The same conclusions can be drawn from Fig. |4|(g) 
and (h) which show the LDOS at the Fermi level of a 
p-InSb surface covered with 1.5 % Cs. Again a regular 
wave pattern not exhibiting simple ring like structures 
is visible and the Fourier transform exhibits a ring 
with additional intensity in the center. The resulting 
fc-values of the ring are compared with ARPES results 
and a parabolic dispersion (Fig. |4][h)), both resulting 



Interestingly, the corrugation of the LDOS pattern 
is distributed rather homogeneously across the sample. 
The corrugation strength, defined by the ratio between 
spatially fluctuating dl /dV -intensity and total dl/dV- 
intcnsity according to 



LDOS 



dl/dVm 



dl/dV^ 



dl/dV^ 



(20) 



with dl /dVYnin and dl / dV-axca.i\ determined as sketched in 
Fig. [5]jd), is of the order of 50 ± 5 % as shown by the 
scaled histogram in Fig. |5jb) and (c). This value is 
much larger than the corrugation strength in the 3DES, 
which is 3 ± 0.5 % as shown in Fig. Eld).*^ Moreover, 
the corrugation can be reproduced by the single particle 
calculation according to eq. 12 (see Fig. [5](c)). Both 
results reflect the much higher probability of multiple 
scattering of an electron with in a 2DES which eventu- 
ally leads to weak localization.'i22l Thus, many different 
scattering paths containing each many scattering events 
contribute to the LDOS leading to more intricate pat- 
terns. The fact, that the probability for back scatte ring 
is enhanced in 2D, results in the stronger corrugation.l^J^ 
Since we map the LDOS and not individual states, we 
cannot determine the localization length as long as it is 
larger than the average distance between localized states 
contributing to an LDOS image being about 65 nm in 
Fig. |4]^g). However, the transition between strong lo- 
calization and weak localization can be determined by 
the change of LDOS corrugation This is demonstrated 
in Fig. [5jd) showing histograms for different voltages 
of the Co/p-InAs(110) system. The low dl/dV peak at 
high voltages corresponds to the Co islands and, thus, 
not to the 2DES. Obviously, the 2DES peak starts at 
dl/dV = nS at low energy, but gets distinct from the 
peak of the Co islands at higher voltage. The corrugation 
strength Cldos drops from about 100 % to about 50 % 
at a sample voltage of 50 meV, where the electrons be- 
come free from local confinement, i.e. strong localization. 
Although we do not have a quantitative understanding 
of this drop, it appears to be a strong indication of the 
transition between strong and weak localization. 
Obviously, disorder potential and LDOS of a 2DES can 
be mapped rather precisely. This is a crucial result, since 
the typically known band parameters,!^ the potential 
landscape, and the electron density completely determine 
the appearance of the LDOS. Thus, detailed comparison 
with theoretical models becomes possible. 
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FIG. 5. (a) histogram of dl/dV values obtained on n-InAs(llO) (A^d ~ 1-1 ■ lO'^^ m ^) after cleavage with corrugation Cldos 
marked {V = 0.05 Y, I — 0.4 nA, V^od = 8.5 mV), inset: part of the dl/dV image used for evaluating the histogram; (b) same 
as (a) for n-InAs(llO) (A^d = 1.1 ■ 10^^ m"^) covered with 2.7 % Fe {V = 0.10 V, 7 = 0.3 nA, V;^od = 1.8 mV); (c) histogram 
of dl/dV values of n-InAs(llO) (Ad = 1.1 • 10^^ m"^) covered with 0.8 % Fe {V = 0.10 V, 7 = 0.3 nA, V^od = 1.8 mV) 
in comparison with the histogram of calculated LDOS by diagonalizing the Hamiltonian with the experimentally determined 
disorder potential using periodic boundary conditions; the experimental curve is stretched by 5 %; (d) histogram of dl/dV 
values obtained on p-InAs(llO) (TVa = 4.6-10^^ m"^) covered with 15 % Co at the voltages marked (K,tab = 0.43 V, 7stab = 0.7 
nA, Knod = 5 mV); mean and min as used for the determina tion o f Cldos are marked; (e) corrugation Cldos as a function 
of energy deduced from the histograms shown in (d); T — 6 K! ^^ * ^" * 



VI. LOCALIZED AND EXTENDED STATES IN 
THE QUANTUM HALL REGIME 

One of the most surprising effects of a 2DES is the 
quantum Hall effect (QHE)-'^, which is often used to 
determine the quality of a 2DES. Applying a magnetic 
field -B perpendicular to the 2DES and driving a current 
/ through the 2DES, one finds certain i?-field regions, 
where the voltage drop parallel to the current vanishes, 
i.e. Vxx = V, while the voltage transversal to the cur- 
rent direction Vfjaii exhibits plateaus exactly at 

l^Hall ^^-I (21) 

where i, the so called filling factor, is the integer clos- 
est to h ■ -/V2DEs/(e7?). The value of Ru — Vnaii/-^ is 
based on fundamental constants only and independent of 
the details of the 2DES. It is used, e.g., as a resistance 



standardiiiJ and di scus sed as an ingredient with respect 
to a mass standard.l^i^ 

The Hall resistance of an adsorbate induced 2DES is 
shown in Fig. |6][a) for different electron densities with 
the plateaus for i = 1 — 4 marked. 

The origin of the QHEisan interplay of Landau quan- 
tization and disorder .lii^llli^ Without disorder the single 
particle energies of the 2DES are quantized according to 

Em^n,s = (n - 0.5) • ^^"j^ + s ■ g{E)tiBB + (22) 

with s = 1/2,-1/2 being the spin quantum number, n, 
m being positive integers and fiB — 5.8 • 10^^ eV/T be- 
ing the Bohr magneton. The first two terms are dubbed 
Landau quantization and spin quantization and Em is 
the subband energy of the 2DES discussed in chapter jll] 
The development of these energies with B field can be 
probed by STS as shown in Fig. [6]jb). The measured 
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FIG. 6. (a) Hall resistance of p-InSb(llO) (Na = 1 — 2 ■ 10^^ m^"^) covered with Ag of the amount indicated, the resistance 
is measured by 4-probe lock-in technique (13 Hz), T = 2 K (courtesy of R. Masutomi, Tokyo (b) grayscale plot of dl/dV 
intensity as a function of voltage and applied magnetic field B, green and yellow lines mark spin-down Landau levels of first 
(El) and second (-E2) subband, respectively, Landau level (LL) numbers and spin directions of the first subband are marked 
on the right (Ktab ~ 0.15 V, /stab ~ 0.1 nA, Vmod = 1-5 mV); (c) grayscale plot of dl /dV intensity as a function of voltage and 
tip position with respect to the surface (B — GT, Ktab = 0.15 V, /stab ~ 0.13 uA, Vmod ~ 1-3 mV); (d)-(j) dl/dV images at 
B = 12 T and at the voltages indicated in (1): V = -116,3 mV, -111.2 mV, -104.4 mV, -100.9 mV, -99.2 mV, -92.4 mV, 
—89.0 mV, i.e. across the lowest spin-down Landau level, white arrows in (b), (e), (i), (j) mark drift states in potential valleys, 
green arrows in (h), (j) mark drift states on potential hills, red and yellow arrows in (f), (h) mark remaining intensity at saddle 
points; crosses in (g) mark saddle points, too (Vstab ~ 0.15 V, /stab = 0.1 nA, V^od = 1-0 mV); (k) calculated LDOS at the 
center of LLO ^ at /? = 12 T, crosses and red arrows have the same meaning as in (f)-(h); (1) spatially averaged dl /dV curve 
corresponding to the images in (d)-(j) as marked, arrows mark the spin direction of the lowest Landau level LLO; (b)-(j),(l): 
r = 0.3 K.''^ 



dl /dV intensity at a single point is shown as a function 
of V and B exhibiting two fan-like ensembles of lines, 
which are labeled by subband energies E„i. The differ- 
ent Landau levels (LLn) and spin levels (t, |) are visible 
and the linear energy dependence on B field is discern- 
able. The corresponding single particle states without 
disorder are highly degenerate and can be described as 
rings with different radii for different LLn encircling a 



single flux quantum = h/e each. This leads to a level 
degeneracy of Nn^m,s = eB/h independent on the Lan- 
dau level number or the spin direction. 
However, these states are subject to random potential 
disorder within a semiconductor. Thus, they change their 
energy as a function of position as demonstrated in Fig. 
|6]jc), where the pairs of Landau levels belonging to dif- 
ferent spin quantum numbers meander in energy as a 



14 



function of position. The fact that the meandering am- 
phtude is larger for the lower lying Landau levels mimics 
the fact that the radius of the rings called cyclotron ra- 
dius Ten = \/ (2(n — 1) + 1) • ?i/ (eB) increases with the 
Landau level index. Thus, the LDOS of higher Landau 
levels probes the potential on a rougher length scale.'^^ 
In 2D, the disorder potential has primarily a semiclas- 
sical effect i^^'^ the electrons perform the fast cyclotron 
rotation within the electrostatic disorder, which leads to 
additional drift moti on a long the equipotential lines of 
the disorder potentiaLl^^ Basically, the electrons are ac- 
celerated and decelerated, if they move downhill or uphill 
within the potential disorder during their cyclotron ro- 
tation. This results in different radii of curvature of the 
electron path at lower and higher potential energy, di- 
rectly implying a motion perpendicular to the gradient 
of the potential as long as the gradient direction remains 
similar on the length scale of the cyclotron radius. Quan- 
tum mechanically, so-called drift states meander along 
equipoten tial line s with a width of about the cyclotron 
radius re„.™lll 

If the potential energy of the state is low (high), the drift 
states are closed trajectories around potential minima 
(maxima), i.e. they are localized and represent insulat- 
ing electron phases. Thus, whenever if the Fermi level is 
located at energies belonging to localized states, the lon- 
gitudinal conductivity a^x vanishes for T — >■ K. Since 
the longitudinal resistance p^x is: 



Pxx 



(23) 



and the Hall conductance is a^y ^ S at i? > T, p^x 
vanishes as well. 

Only, in the center of a LL, the equipotential lines of adja- 
cent potential valleys and hills merge at the saddle points 
of the potential leading to an extended state. It can be 
shown that exactly one equipotential line travers es th e 
whole sample in the limit of infinitely large samples.'^' It 
is known that the state meandering along this line is the 
quantum critical state of the integer QH transitions and 
responsible for the finite longitudinal resista nce be tween 
quantized values of the Hall conductance^^'' H^'^i^. This 
quantum phase transition between localized states in the 
valleys via an extended state towards localized states at 
hills of the potential disorder is repeated for each (m, n, s) 
level of eq. 
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This transition is universal, i.e. the 
energy dependent localization length £,{E) of the states, 
defined as 

< \'i>^{x - x^^d?) ■ 5{E ~ E,) >,oc e-l-l/«(^) (24) 

{xq j! center of mass of l^'ip) is independent of details of 
the disorder being 



^{E) ^\E- E, 



(25) 



Vc is not known analytically, but has be en evalu ated nu- 
merically, which results in Vc — 2.4— 2. 6. ^^^^^^^ Another 
univ ersality is the multifractal spectrum of the critical 
state.™Sl 

The corresponding transition is shown for the LDOS of 
the lowest energv (m, n, s) level of an adsorbate induced 
2DES in Fig. [6fd)-(l).l^ The spatially averaged dl/dV 
curve is shown in Fig. |6][1) exhibiting the single peak of 
LLI4,. The energies of the LDOS images are marked. In 
the low-energy tail of the peak (Fig. ^&)), the LDOS 
exhibits spatially isolated closed-loop patterns with av- 
eraged full width at half maximum (FWHM) ~ 6.9 nm 
close to the cyclotron radius — 7.4 nm. These are 
localized drift states of the n ~ l,m = 1,S = —1/2 
level LLO4, aligning along equipotential lines around a 
potential minimum. Accordingly, at slightly higher en- 
ergy [Fig. [6](b)], the area encircled by the drift states 
increases indicating that each drift state probes a longer 
equipotential line at higher energy within the same val- 
ley. In contrast, the ring patterns at the high-energy 
tail, marked by green arrows in Fig. ^i), (j), encircle 
an area decreasing in size with increasing voltage. These 
states are attributed to localized drift states around po- 
tential maxima. Notice that the structures in Fig. l6td) 
and (e) appear nearly identical in Fig. and (j) as 

marked by white arrows. The latter structures are the 
LLlf states localized around potential minima, which en- 
ergetically overlap with the high-energy LL1| states lo- 
calized around potential maxima. When the voltage is 
close to the LLlj, center [Fig. |6]Jf), (h)], adjacent drift 
states coalesce and a dense network is observed directly 
at the LL center [Fig.[6](g)]. This is exactly the expected 
behavior of an extend ed drift state at the QH transi- 
tion as described above.l^^^E^Fig. [6|^k) shows the calcu- 
lated LDOS around an extended state at J5 = 12 T in 
a 2DES of InSb with a random distribution of dopants 
with the experimentally known densities = 9 • 10^^ 
m~^ and = 5- 10^^ m"''. Each dopant as usual is sur- 
rounded by a scre ened Coulomb potential. The Hartree 

125 



with E„it being the energy of the extended state and 
i/c being the universal, critical exponent. The value of 



approximation ■ ' is used for the diagonalization of the 
matrix. Good qualitative agreement with the measure- 
ment is achieved supporting the interpretation of the co- 
alesced LDOS patterns as caused by the presence of an 
extended state, albeit the disorder in the calculations is 
obviously smoother, which could be traced back to a re- 
maining influence of the Cs chains thatpartially form on 
the surface.Esl The crosses within Fig. ^g) and (k) mark 
the saddle points of the potential appearing as extended 
LDOS areas due to the vanishing steepness of the poten- 
tial at the saddle points. The same points are marked 
by arrows in Figj6|^f) and (h), where remaining intensity 
is found at the saddles. It has been argued initially that 
this intensity is related to the quantum tunneling across 
the saddles w hich is decisive for the value of the criti- 
cal exponent t/l ii5 | i22 |^^.j . im nierical calculations revealed 
that this is not the case^^ Instead, it reflects only the 
energy resolution of the experiment of 2.5 meV being sig- 
nificantly larger than the expected SE ~ 3.3 • /cbT = 0.1 
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meV.'fifi'This might be related to the slow Coulomb glass 
dynamics in localized systems p^ ^ ^^'' ^ but further experi- 
ments are required to clarify this issue. 
Similar transitions between localized and extended states 
have also been observed on the surface of highly ori- 
ented pyrolytic graphite (HOPG^pS! and on turbostr atic 
graphene on top of graphite on SiC(000T)p2HllI] Re- 
cently, it was also investigated in higher Landau levels 
of the InSb 2DES, where an additional nodal structure 
of the LDOS perpendic ular to the equipotential lines was 
revealed and explained. ^'^^ 

Importantly, the LDOS across a universal quantum phase 
transition, the Quantum Hall transition, can be directly 
mapped within the adsorbate induced 2DES. Improved 
theoretical tools using vortex states within the framework 
of real-time Green's-functions allow, on the other hand, 
a direct calculat ion of LD OS properties for sufficiently 
smooth disorder.fii^E^llIlll Thus, a one-to-one compari- 
son of LDOS patterns pinpointing to the local signatures 
of the influence of particular interactions appears possi- 
ble. 

Of course, the quantum phase transition explains only 
the vanishing conductivity a^x at most B fields but not 
the quantization of the Hall conductance. This can ei- 
ther be explain ed rather elegantly in terms of Chern 
numbera^SMlM] qj. l ess a bstract by the presence of the 
so-called edge states.^^ The latter explanation is based 
on the fact that the drift states at the edge of the sam- 
ple arc meandering along the edge and, thereby, connect 
different leads. Moreover, they are chiral, i.e. they are 
allowing electron motion only in one direction and not in 
the opposite one, which corresponds to the drift direction 
explained above. If the Fermi level is within the local- 
ized states of the Landau level, the edge states at the 
Fermi level are the only ones which can carry the cur- 
rent. Since backscattering is not possible within these 
states, no voltage drop can appear along the edge states. 
Moreover, only a voltage imbalance between states on the 
left hand side of the sample carrying the forward elec- 
trons and on the right hand side of the sample carrying 
the backward electrons can lead to an effective current 
I through the sample. It is straightforward to calculate 
that the relation between effective current / and voltage 
imbalance Vnaii is //VHaii = /h for each pair of edge 
states. Moreover, each critical state below the Fermi level 
requires one of these edge states crossing the Fermi level 
at each side due to the rising confinement potential at 
the edge. These two facts together lead to 
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Hall 



(26) 



with i being the number of critical states below E-p. 
So far the edge states have not been probed by STS, 
but by other scann ing p robe techniques as scanning gate 
microscopy (SGM) j^^^ electrostatic force micr oscopy,^^ 
scanning capacitance microscopy ( SCM ) scanning 
near-field optical microscopy (SN OM) and a scanning 
single electron transistor (SSET).'i3Il However, all these 



techniques exhibit a much worse lateral resolution well 
above 100 nm. Thus, it would be interesting to probe 
the edge states by STS, too. SSET has also been used 
to map properties within the fractional quantum Hall 
regime, which is guided by electron-electron interactions. 
These experiments showed rat her directly the fractional 
excitation charge of AQ = e/S.'SH 



VII. PROBING THE EXCHANGE 
INTERACTION 

One of the most simple type of electron-electron inter- 
action is the exchange interaction. For a system of two 
particles, it reads 

(27) 

where Vcd^x) is the electron-electron interaction poten- 
tial. 

Due to symmetry of the many particle wave function, 
the exchange interaction is attractive for parallel spins 
and repulsive for antiparallel spins. Thus, an electron 
within a spin-polarized 2DES exhibits a lower effective 
electron-electron repulsion energy, if parallel to the spin- 
polarization, than an electron of the same spin that is 
embedded into a spin neutral 2DES. An electron with 
antiparallel spin to a spin-polarized 2DES exhibits an 
even stronger electron-electron repulsion. 
Applying a magnetic field to a 2DES leads to an os- 
cillating spin polarization. In case, the filling factor 
V — hN/{eB) is odd, the spin polarization is maximum, 
i.e. only the half of the highest Landau level with the 
spin moments antiparallel to -B is filled with electrons. 
This leads to a spin density Ng = s- eB/h or to a magne- 
tization M ~ hb ■ eB /fi-b/i = 1.4 A/m-B[Tesla]. In case 
of an even there is no spin polarization, since the same 
amount of f levels and \. levels are filled with electrons. 
In the case of spin polarization, the lower energy spin 
levels of a Landau level will gain energy with respect to 
the ones in the higher energy spin level. Thus, the split- 
ting of the two spin levels will be larger than g{E)fiBB 
in the case of spin polarization, and exactly g{E)fiBB 
without spin polarization, i.e. at even v. This leads 
to an oscillatory effective ^eff (^)-factor determining the 
spin splitting Ai?ss(i3) = QesiB) ■ hb ■ B as a function 
of B. The increased spin sp litting at odd v is called ex- 
change enhancement lUSHifil 

Within the QH regime, exchange enhancement is a rather 
local effect, since the expression in eq. [27] requires an 
overlap of different localized wave functions. Conse- 
quently, the exchange enhancement depends on the lo- 
cal spin polarization mostly within a single valley or hill 
region of the potential disorder. Exchange enhancement 
has, however, previously only been measured without lat- 
eral resolution, e.g. by capacitance spectroscopy-fSS 
Figure T^a) shows the band structure of an adsorbate 
induced 2DES on p-InSb(llO) at low doping. One ob- 
serves that the band bending reaches about 600 nm into 
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FIG. 7. (a) Band bending as calculated by a Poisson-Schrodinger solver for p-InSb(llO) (A^a 
% Cs, valence band maximum (blue), conduction band minimum (red), Fermi level -Ef and |^n,(z)p of the first two subbands 
are shown in the direction perpendicular to the surface, inset: magnification of the near surface area; (b) grayscale plot of 
d//dV-intensity as a function of magnetic field B and applied sample voltage V; Fermi level (dashed line) and two spin levels 
of LLO are marked {T — 5 K, 14tab = 0.3 V, /stab = 0.4 nA, Knod = 1-6 mV); (c) energy splitting between t and ^ level of 
the lowest Landau level LLO (black line) in comparison with the splitting expected at \g(E)\ = 42 (red line), inset: dl/dV 
curve (black) in the region of LLO in comparison with two Gaussians colored yellow and leading to the transparent red dl/dV 
curve, B — 5.5 T; (d) difference between black and red curve from (c) (black line), smoothed curve is shown in green, scale 
bars mark the exchange enhancement at the corresponding filling factors calculated within the random phase approximation; 
(e) colorscale plots of dl/dV intensity of graphene on SiC(OOOl) as a function of voltage and tip position with respect to the 
surface for different magnetic fields B and local filling factors v as indicated, levels corresponding to different spins t and J, and 

K/K' combin ations are marked, (T = 0.01 K, Vstab = 0.25 V, /stab = 0.2 nA, Vmod = 0.05 mV) (courtesy of J. Stroscio, NIST 
Gaithersburg).E2Iill 



the bulk of the sample. This implies an insulating region 
of more than 500 nm thickness between the 2DES at the 
front of the sample and the conducting bulk region of the 
InSb. Consequently, the screening of the bulk electrons 
can be neglected on the scale of the electron-electron dis- 
tance within the 2DES being only about 10 nm. This 
increases the strength of Vcc- Additionally, the insulat- 
ing barrier cannot be penetrated by electrons, thus the 
current between tip and 2DES must flow through the 
2DES, which had to be contacted on the side by a silver 
wire.l^ Since localization reduces the conductivity of the 
2DES in the QH regime, one has to be careful that the 
tunneling electrons, which tunnel from the 2DES to the 
tip, are replaced in between individual tunneling events, 
i.e. the current must be low enough such that local equi- 
librium is probed by each tunneling event. Otherwise, 
the remaining local charge will give rise to an additional 
local band bending and, thus, to a locally modified and 
time-dependent filling factor. Experimentally, we found 
that a tunneling resistance above 1 Gft could fulfill this 
requirement up to i? = 7 T and down to T = 5 K, i.e. 
the band bending by this so-called spreading resistance 
remains significantly less than the distance between ad- 
jacent levels.*^ 

Figure ^h) shows dI/dV{V) spectra of such a 2DES 



taken at a fixed position while ramping the magnetic 
field B. A Landau fan similar to Fig. [sjb) is visible. 
However, in contrast to Fig. [sjb), the lines of conduc- 
tance maxima are wavy and not straight. The reason for 
waviness is the fixed electron density A^2DES which has to 
adapt to the degeneracy of the Landau levels. Thus, the 
distance between Ep and the lowest Landau level jumps 
each time, Ep has to move into the next lower Landau 
level. Since Ep is fixed in the tunneling experiment a 
wavy movement of all the Landau levels with respect to 
Ep results. 

More importantly, the distance between adjacent spin 
levels AEss deviates from g{E)fiBB. To analyse this 
effect, barely visible in Fig. [7]Jb), in more detail, we con- 
centrate on the lowest LL around —120 mV, which gives 
the highest accuracy in determining Ai?ss- We adapted 
two Gaussians for all 386 spectra between 3.5 T and 
6.1 T. The Gaussians, having equal width and height, 
were fitted using a nonlinear least squares method and a 
trust-region algorithm as implemented in Matlab.^"*^ The 
fits are good as can be seen in the inset of Fig. ^c) and 
by the large confidence value of = 0.94 (0.97 above 
5 T). The error for the resulting AEss is about 0.2 meV. 
The resulting spin splitting AEss{B) as function of B is 
shown in Fig. u\c) in comparison to a straight line cor- 
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responding to ordinary Zeeman splitting of \g{E)\^BB 
with \g{E)\ = 42. 

Figure [7|[d) , finally, shows the deviation A{B) from the 
straight line. It oscillates around mcV with maxima 
(minima) around odd (even) filling factors as expected 
for the exchange enhancement. Negative values of A(i?) 
are probably caused by slight deviations from a spin split- 
ting linear in B due to increased spreading resistance 
with increasing i?, which leads to superlinearity, and non- 
parabolicity of InSb leading to a smooth decrease of g{E) 
with B, thus, supralinearity. However, both effects are, 
in first order, monotonic in B and cannot explain the 
oscillations. One could imagine that the spreading re- 
sistance depends oscillatory on filling factor. But then, 
it would be largest at even filling factors, where Ep has 
the largest energy distance to the next critical state. This 
would lead to an oscillation of A{B) with maxima at even 
filling factor in contrast to the experimental observation. 
Thus, the oscillatory B field dependence of the spin split- 
ting is in accordance with an exchange enhancement of 
about 0.5-0.7 meV. 

To substantiate this assignment, we calculated the ex- 
pected exchange enhancement in the lowest Landau level 
using a random phase approximation (RPA). This ap- 
proximation neglects the dielec tric scre ening at other fre- 
quencies than the exciting one.'— -'--P' This is well justi- 
fied since the electron density in the 2DES Af2DES is large 
compared to the scale set by the eff ective Bohr radius, 
i-e. iVzDES - 25- (m*/(eaB) respectively, the r^- 
parameter is much smaller than one. We performed the 
calculation using m* = 0.02 and g{E) = —42 as deduced 
from the Landau level distance and spin level distance 
visible in Fig. [7](b). 

The 2D bare Coulomb potential represented in Fourier 



space is given by 



V{q) = F{q) 



eq 



(28) 



with the form factor F{q) accou ntin g for the finite exten- 
sion of the 2DES in z directionpsl 

^(^) = £ + 8i + 4^-^*^ (2^) 
X = 1 + ^ . (30) 

3/ 4S7Tm*mce^ (11 



(M-^2DES + ^Dopant) 



with the dielectric constant of InSb e — 16.8, and 
iVzDES = 2.7 X 10^6 m-2 as well as iVoopant = 8 x 
10^"* m^^, the 2DES density and the density of ionized 
acceptors, respectively. 

Restricting ourselves to the static response of the 2DES 
to perturbations, the dielectric screening of the potential 
requires to replace V(q) by 



VRPA(g) 



V{q) 



with 



l~V{q)Ii^{q) ' 



(31) 



(32) 

Here, f{E) = (1 + pJ^^^^j/lfeBT))-! ^re Fermi functions 
evaluated at T = 5 K, En,s are the energies corresponding 
to the n. Landau level and spin level s, and 

Pn,n.{x) = i-ir+"^e-^U:,-"^x)L^-^\x) , (33) 

with L™^"(a:;) being associated Laguerre polynomials. 
The total exchange energy for an electron in Landau level 
n with spin s then reads 



qdq 
2tt 



^ till - Em,s + 



(34) 



where fi; = (21 + 1) -n ■ {k-gT) (1: integer) are Matsubara 
frequencies. 

The filling factor dependent exchange enhancement 
EEn{v) within the first Landau level [n = 1) is then 



(35) 

This leads numerically to, e.g., EEi{v — 11) = 0.76 meV 
and EEi{i' = 15) = 0.55 meV. The calculated EEi[v) 
for V = 15,13,11,9 are added as vertical bars in Fig. 
[Tjjd) showing very good agreement with the experimen- 
tal data. It has been checked that the static screening 
approximation used here is suffici ent to reproduce the 
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results of more elaborate approache*'^ on a quantitative 



level. Moreover, the results for EE^) barely change, if 
the system parameters g{E) and m* (respectively N) 
are varied within reasonable limits. For example, the 
change of EEq is less than 1 % for g{E) = —38 and still 
below 10 %, if m* is doubled to m* = 0.04.1^3 
The fact that magnitude and oscillation phase of A(i?) 
compare favorably with a parameter free calculation 
of the exchange enhancement, is strong evidence that 
this short-ranged electron-electron interaction effect is 
quantitatively probed by STS, despite the screening 
properties of the close-by tip. Notice that the potential 
disorder and the second subband of the 2DES are 
not included within the calculations and, thus, are of 
minor importance for the exchange interaction which 
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happens on the length scale Vc — 10 nm smaller than the 
correlation length of the disorder potential (30-50 nm). 
The same 2DES also exhibited a Coulomb gap at Ep, 
which is a sign for the presence of long-ranged Coulomb 
repulsion. ^^^^^'^ Here, however, the screening by the 
tip had to be taken into account, in order to describe 
the observed Coulomb gap quantitatively,!^ i.e., the 
tip is described as a screening gate 8.6 nm away from 
the plane of the 2DES. This stresses the need for 
less screening tips, as e.g. tips made from low-doped 
semiconductors 

Exchange enhancement has also been obse rved on the 
system graphene on graphite on SiC(0001).'i^ At very 
low temperature (T = 13 mK), graphene shows spin 
splitting and valley splitting. This is demonstrated by 
the color plots in Fig. [7]je), where four lines belonging 
to a single Landau level meander in energy as a function 
of position. Notice that the plots are rotated by 90° 
with respect to Fig. ^c), i.e. position is along the 
vertical and energy along the horizontal axis. Moreover, 
plots recorded at different B, as marked on the right, 
are plotted above each other. At i? = 12 T, one pair 
of lines separated by the Zeeman energy of 1.4 meV is 
observed above and one pair of lines appears below 
the Fermi level. Each pair corresponds to one of the 
sublattices of graphene at the measurement position, 
i.e. to a different, mutually orthogonal, combination 
of the two valleys K and K' as marked. All four lines 
belong to the n = 1 Landau level. By decreasing the 
magnetic field, the degeneracy of each line ull decreases 
according to ull = eB/h and the levels, which are above 
Ep at higher B must cross E-p. At the crossing point of 
the first spin level, the distance between the spin levels 
increases by about 4 meV and it decreases again, if both 
spin levels have crossed Ep. This is an even clearer 
example of the exchange enhancement, which profits 
from the low bare g-factor g — 2 of graphene making the 
bare spin splitting smaller and of the smaller extension 
of the graphene electrons in z direction enhancing their 



and eq. 29 



Careful inspection of Fig. [vje), reveals a 
weak third line (marked by arrows) at half-valued filling 
factor, which is interpreted as an interaction effect with 
the quan tized 2D system of the underlying graphene 
layer.^ 

The experiments on exchange enhancement so far have 
not been used to probe local differences of the exchange 
enhancement at the same B, which is the obvious 
strength of STM, and, thus, outlines an interesting 
subject for future research. A prototype experiment has 
been done using the spin splitting of states within the 
TIQD, where a non-local correlation of the exchange 
enhancement with the disorder potential has been 
found.lail 

VIII. PROBING THE LOCAL RASHBA EFFECT 

As mentioned in the introduction, the exchange in- 
teraction between electrons is used within spin based 
qubits,'^ which allows to manipulate the spin degree 
of freedom electrically via the tunable overlap of wave 
functions in adjacent quantum dots. For propagating 
electrons, also the spin-orbit coupling, which eventually 
is a relativistic effect originating from the Dirac equa- 
tion, can be used to manipulate spins via electric fields.l^ 
The Rashba term within the Schrodinger equation can 
be written asP^ 



Rashba 



aE ■ {kx a) 



(36) 



Coulomb interaction significantly according to eq. 28 



with k being the wave vector of the electron, a being 
the Pauli matrices and E being the unit vector along the 
electric field typically perpendicular to a surface or an 
interface. In first order, the Rashba-parameter a within 
III-V semiconductors depends on the strength of the elec- 
tric field E, the energy gap i?Gap) the effective mass m* 
and the spin-orbit splitting Aso within the valence band, 
which is large for crystals made of heavy atoms. A good 
approximation b ased on an 8-band k-p description of the 
band structure isjlSlIIIZl 
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Since In and Sb are heavy 5p-atoms, they exhibit a large 
Aso = 0.8 eV at T = 4 K. Moreover, as described in 
chapter [TT| and [V| band gap Eq^p = 0.23 eV and effective 
mass at the band edge mg = 0.0135 are quite low, such 
that a relatively large prefactor of a with respect to E 
results: 



5.1-10-18 



eVm^ 
V 



\E\ 



(38) 



Notice that the spin-orbit term Aso within binary, 
ternary or quaternary alloys of Bi, Sb, Se, and Te can 



even lead to a band inversion, i.e. part of the p-type 
band becomes the conduction band while the part of the 
s-type band becomes a valence band. Thi s leads to the 
currently celebrated topological insulators liSMlSlI -^vhich 
exhibit, at least, one surface state with spin chirality 
within the fundamental band gap, which, moreover, is 
protected by time-reversal symmetry. 
In order to get a strong a, one needs, in addition to the 
prefactor, a strong electric field. An interesting ques- 
tion is, if charge neutrality prohibits the presence of 
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1.5% Cs/p-lnSb (110), Na= I.S IO^" m-^ 
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FIG. 8. (a) Band bending as calculated by a Poisson-Schrodinger solver for p-InSb(llO) (TVa ~ 1.5 ■ lO'^* m~^) covered with 
1.5 % Cs, conduction band, valence band, and band gap are marked, of the first two subbands are shown in yellow 

and electric field-B is marked in red, field value and resulting Rashba parameter a are given; (b) dl/dV curves recorded at 
different positions of the sample, spin directions are marked (_B = 7 T, I4tab = 0.3 V, /stab = 0.2 uA, Knod = 1-0 mV); (c) 
relation between the peak voltage of the lower spin level and the determined spin splitting; rings mark two different areas 
where spin splitting is observed and indicate that it is only observed in hill and valley regions of the potential; (d)-(i) spatially 
averaged dl/dV curves (area for averaging: (300 nm)^) at S-fields indicated (Ktab = 0.3 V, /stab = 0.2 nA, V^od = 1.0 mV) 
in comparison with calculated DOS using m* — 0.035, a = 7 ■ 10~^^ eVm, g{E) — —21, and a Gaussian level broadening with 
FWHM as marked, double arrows mark positions of nodes in DOS and dl/dV curve; T — 5 kW^ 



electric fields within a 2DES. Winkler has shown that 
the relevant property for a is the effective field within 
the valence band, which remains present due to different 
band parameters in valence and conduction band, and 
that this more correct description is except for prefac- 
tors close to on e ide ntical to the quantitative descrip- 
tion given ahave^^ The strongest field is found at the 
surfaces of metals (10^ — 10^ V/m). Consequently, the 
strongest a have been found for surface states on met- 
als, e.g.. for Bi(lll), one has found a = 0.55 • lO"^'^ 
eVni,liS^ while for Bi alloys even a — 3.05-10"^° eVm has 
been obtained.l^^ The latter value is partly attributed to 
additional interatomic electric fields. ^^'^ However, within 
metals the electric field can barely be controlled by gates. 
Thus, spins can not be manipulated on ns time scales as 
required for information processing. Most likely, InSb is 



the semiconductor with the largest possible prefactor of 
a not being a topological insulator. 

The first consequence of Hso is a modification of the 
band structure. Assuming a parabolic dispersion of the 
kinetic energy as for the conduction band of InAs and 
InSb, the origin of the parabola in k is offset from the 
F-point, which marks the center of the Brillouin zone, by 

|Afc| = (39) 

The fc-direction of offset is different for the states with 
spins pointing to the right with respect to k and for the 
spins pointing to the left with respect to k. Thus, spin 
degeneracy is lifted and the eigenstates of the spin are 
always oriented perpendicular to k and The resulting 
two parabolas for one k direction are rotationally sym- 
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metric leading to a circular trough in the E(k) dispersion 
with radius | Afc|. A cut through k space in E{k) using the 
plane perpendicular to E and energies above the bottom 
of the trough exhibits two circles with chiral spin texture 
perpendicular to k. Opposite spin directions distinguish 
the two different circles. The difference between the radii 
of the circles is exactly the |Afc| of eq. 39 i.e. indepen- 
dent of energy as long as m* is independent of energy. 
The Rashba-split E{k) dispersion has been m easured 
directly by ARPES, but so far only for metalP^HUl] 
or for metal/semiconductor interface states f^SSHlZSI ]^■^^■^ 
not for intrinsic states of a semiconductor. The rea- 
son is the low m* , e.g. for InSb, that results in small 
|AA:| ~ 1.3/V • \E\ between the two circles. This requires 
excellent angular resolution in ARPES of the order of 
0.1°. Instead, the Rashba effect in Ill-V-semiconductors 
has been probed b y the be ating pattern of Shubnikov- 
de-Haas-oscillation^^SEZll or by the analysis of weak 
antilocalization.^^S'i^ The beating can be explained by 
the two different circles in k space at the Fermi level, 
which are both not spin degenerate. Since each cir- 
cle with radius |fc|F,i (* = lj2) encircles a density of 
= l^lFi/(4^) ^-points and the degeneracy of a spin 
polarized Landau level is eB/h, the number of filled Lan- 
dau levels (within the circle) is: 



AireB 



(40) 



Each time, when i^i is an integer, the Landau level high- 
est in energy is completely filled and the conductivity 
exhibits a minimum. In contrast, a half-valued Vi, i.e. a 
Fermi level in the center of a Landau level leads to a max- 
imum in conductivity. If, both, i>i and i>2 are integers the 
corresponding minimum is deep, while if vi is an integer 
and i/2 is half-valued, the minimum disappears. From 
this beating of the oscillation in conductivity one can de- 
termine |Afc| and, thus, 

The same effect leads to a beating of the density of states 
which can be probed by STS. A node of the beating will 
be obtained at an average |A:|-value of the two circles A:„, 

respectively an energy i?„ — fi^ k^ / {2mem*) with respect 
to the subband energy according to: 



n + 0.5 = 1^1 - V2 (41) 
^ h(kn + |Afc|/2)^ - (fc„ - |Afc|/2)^ 
AneB 

2hkr,\Ak\ 
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eB 

h\Ak\ ' 

h*e^B^ 
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heB 



(n-HO.5) 



(42) 
(43) 
{n + 0.5)(44) 

(45) 



Using a = 10"^° cVm, S = 6 T, and m* = 0.03, one 
gets En = 150 meV -(n -|- 0.5), i.e. the first node of the 
beating must be about 75 meV above the subband onset. 



Notice that the energy position of the beating node de- 
pends quadratically on i3-field and quadratically on the 
inverse of the Rashba parameter a. 

A more sophisticated description considering Zecman 
splitting and Rashba spin splitting on equal footing has 
been given by Rashba in his first publication on the 
effect P 



e: 



Ej + TiLo, 



1/2 



^ — a {2m* /fC'uj^ ^ , 
^^lL_m*giE)\ 
2 V 2m5 ; 



(46) 
(47) 
(48) 



Here, n = 0, 1, 2, . . . is the Landau level index, s = 1/2 
for n — 0, and s — ±1/2 for n = 1, 2, 3, . . . being the spin 
index, i is the subband index and lOc — eB / (m* rric) . 
Fig. |8ja) shows the Poisson-Schrodinger result for a 
2DES, which is induced by 1.5% Cs on the strongly 
p-doped InSb(llO). The strong doping results in a 
strong electric field (> lO^V/m) within the 2DES and, 
thus, according to eq. [38] to a large Rashba coefficient of 
a ~ 10~^°eVm. A more sophisticated calculation taking 
the curvatur e of th e bands and its overlap with '^i{z) 
into account ,1^^^^^ but still neglecting the penetration 
of the ^i(z) into vacuum, leads to a ~ 9 • 10~^^eVm. 
This value is larger than the a-values observed in 
InAs inversion-layers or heteros truc tures by transport 
measurements (3-4 • 10~^^ eV m).'^^ It should be noted 
that the calculated Rashba parameter is an upper 
estimate because of the ignored barrier penetration of 
the electronic wave functions. Furthermore, a is only 
the lowest order of an inversion asymmetry induced spin 
splitting and it is known that higher orders lead to a 
reduced effect .1^^ 

Figure [8][d)-(i) shows spatially averaged dI/dV{V) 
spectra of the Cs covered InSb(llO) at differnt B. The 
curves are averaged from 144 curves recorded on a 
regular grid covering an area of (300 nm)^. The onset of 
the first subband at about —50 mV and the onset of the 
second subband at about 150 mV can be identified as 
steps in the dl / dV signal. The energies are in reasonable 
agreement with the result from the Poisson-Schrodinger 
equation. Obviously, only the first subband is occupied 
by electrons. The electron concentration of the 2DES is 
A^2DES — 6.5 X 10^^ m~^. On top of the steps oscillations 
are visible due to Landau quantization, which feature, 
in addition, a changing amplitude. The minima of the 
oscillation are marked by arrows. The distance of the 
peaks is in agreement with the expected distance of 
Landau levels. Thus, the spin splitting is not visible in 
the spatially averaged curves. This is different for single 
curves as shown for two examples in Fig. [sj^b), which 
also shows clearly that the energy shift by disorder of 
about 20 meV is too large to observe spin splitting in 
the averaged curve. Importantly, the beating of the 
spatially averaged Landau level intensity is quantita- 
tively reproduced by the calculations described in eq. [46| 
using a = 7 ■ 10"" eVm, m* = 0.035, g{E) = -21 and 
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a disorder broadening of the levels as marked in Fig. 
|8]jd)-(i). Thereby, m* and g{E) have been calculated 
as an average within the first subband up to the onset 
of the second subband from eq. [17] and eq. [18] The 
parameter FWHM is in excellent agreement with the 
strength of the disorder potential being about 25 meV 
peak-to-peak and shows the expected trend that it 
is coarsened by rd oc 1/\/B. The remaining fitting 
parameter a is only slightly lower than the estimated 
one probably because of neglected barrier penetration 
and higher orders in k ■ p-description as mentioned 
above. The nice correspondence between the measured 
and the calculated node positions at different B is 
strong evidence that the Rashba spin splitting has been 
detected by STM.lSl 

Notice that the Rashba spin splitting is not directly 
visible in quasiparticle interference patterns probed 
by STS, as one would naively expect from the two 
Fermi circles. If only single scattering is considered 
theoretically, the scattering leads to a wave vector 
of the standing waves of exac tly fc, the average of 
the radii of the two circles.fi^ Subtle changes of the 
quasiparticle interference appear, if multiple scattering 
becomes relevant, e.g. within quantum corrals,!^^ and 
a very complex spin distribution within the stan ding^ 
wave should appear, if magnetic scatterers are used.^^° 
Indirect ways to probe the Rashba effect by STS used, 
e.g., the increased density of states at the band ons et of 
the Bi/Ag(lll) and Pb/Ag(lll) surface alloy s.^^U This 
method reveals the Rashba parameter from the strength 
of the peak at the band onset, but it might be very 
sensitive to details of the tip density of states. Another 
method used the complex band structure of Bi(lll), 
which exhibits 6 spin-split Fermi circles, and detects that 
quasiparticle interference requires spin conservation and, 
thus suppresses some of the naively ex pect ed features in 
the Fourier transform of STS images.n^Si This method 
allows to prove spin textures in complex band st ructu res 
and has also been used for topological insulators ,1^^ but 
does not give access to the spin-orbit parameters. 
Local differences of the Rashba spin splitting have not 
been reported so far, althoug h they are of large relevance 
for spin relaxation processes .I^^EllHIHll ^ig. Mb) shows 
that the spin splitting indeed varies with position. For a 
number of curves the spin splitting has been evaluated 
by the same type of fits as shown in the inset of Fig. [tJ^c) 
and the resulting spin splitting is shown as a function 
of the first peak position in Fig. ^^c). Obviously, the 
average spin splitting is larger at higher potential ener- 
gies, i.e. at hills of the disorder potential. The difference 
in average values of the two circles is about 2 meV, 
although the spread is much larger than the uncertainty 
of the individual fits and, thus, depends on unknown 
details of the potential landscape. The stronger spin 
splitting at higher energy has to be contrasted by the 
reduced g{E) with increasing energy. Using g{E) = —31 
as calculated for the subband onset from eq. [17] and 
eq. 18 in combination with the result from the Poisson- 



Schrodinger equation and the average vaules of the spin 
splitting within the left circle at around -57 mV, which 
is AE'gs = 13.5 meV, and the right circle at around -43 
meV, which is Ai?ss = 15.6 meV, we obtain Rashba 
coefficients of a = 5 • 10~^^ eVm and a = 9 • 10~^^ eVm, 
respectively. The extrema of the spin splitting in Fig. 
^^c) require a values of a = 3 — 11 • 10""'^"'^ eVm. Thus, 
the Rashba parameter fluctuates by more than ±50 % 
within the potential disorder. A detailed understanding 
of this fluctuation requires further experiments, but a 
simple estimate offers a tentative explanation for the 
general effect. First, one has to understand that a is a 
local parameter, which depends on the local electric field 
onlyESl Since the band shift at the surface is constant at 
high enough adsorbate density and the z-extension of 

the band bending is proportional to ^ N^l^^^^ (A^2DES 
is negligible at such large doping), the electric field and, 
thus, a is roughly proportional to -/Voopant • Thus the 
question is, how many dopants contribute to the local 
band bending within the area of a single localized state. 
The length scale of the potential fluctuation is 20 nm, 
which corresponds to the lateral extension of localized 
states, while the extension of the electric field in z 
direction is about 20 nm, too, leading to about iVo — 10 
dopants in the electric field area of a single state. The 
standard deviation of zt^/Nu = 3.2 leads to a fluctuation 
of a oc 3.2 ± 0.5, i.e. a FWHM of the a distribution of 
30 %, which is only slightly lower than the experimental 
value of ±25 %. Importantly, already a moderate doping 
leads to a strong fluctuation of a{x,y), which will lead 
to spin dephasing in spintronic devices in depe ndent of 
the well-known Dyakanov-Perel mechanism.lSll 



IX. CONCLUSIONS 

In conclusion, the adsorbate-induced 2DES on 
InAs(llO) and InSb(llO) have been used to detect the 
relevant spin parameters for semiconductor spintronics 
on the local scale. In particular, the exchange splitting 
of electrons has been measured firstly. It depends on the 
local filling factor and is as large as 0.7 meV already at 
rather high filling factors above = 10 for a state far 
away from the Fermi level. Future experiments will in- 
vestigate the spatial dependence in more detail using sys- 
tems with lower filling factor, which naturally would lead 
also to lower distances to Ep and, thus, larger exchange 
energies. Using the conducting bulk of Fig. [7]ja) as an ad- 
ditional gate, would be the best way to reduce the filling 
factor maybe even down to v = 1 ov u = 2. Secondly, the 
Rashba parameter has been measured showing a rather 
large value of a = 7 • 10^^^ eVm on p-doped InSb(llO) 
covered with Cs. Large fluctuations up to ±50 % of the 
Rashba parameter have been observed within a disorder 
potential fluctuating by about 20 meV. This is relevant 
for spintronic devices using the Rashba effect, since it lim- 
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its the effective spin life-time implying a need for more 
detailed studies. Importantly all the STS measurements 
can be done on a 2DES where the potential disorder can 
be measured independently using the tip-induced quan- 
tum dot. Also the other parameters like surface band 
shift and 2DES density can be probed by ARPES. Thus, 
full control on the relevant parameters allows to tackle 
the influence of disorder in great detail. 
Another interesting field with respect to semiconductor 
spintronics are the ferromagnetic semiconductors, where 
room-t empera ture ferromagnetism is still a matter of 
debate.li^MIlSl Mostly, the III-V semiconductors doped 
with Mn have been studied so far by STS, showing, e.g., 
the anisotropic shape of th e Mn acceptor wave function 
resembling a butterfly,U^ its persistence up to nearest 
neighbor Mn-Mn distances, where the local deviation 
from a si mple overlap of butterflies is everywhere less 
than 30 the additional mirror asymmetry of the 

wave function shape appeari ng at d istance up to 8 mono- 
layers from the (110) surfacej ^^'^ l ^^'^ lthe exch ange interac- 
tion between neighboring Mn on the surface,'^^ and the 
critical properties of the wave functions across the en- 
ergy depend ent metal-insulator transition at a Mn dop- 
ing of 1.5 Here, the interplay of disorder and spin 
properties is again very relevant, most likely also for the 



ferromagnetic transition temperatures.'^^ This calls for 
detailed studies by STS also on other types of possible 
ferromagnetic materials with low conductivity. 
Thus, STS probing the spin properties of semiconductors 
will sincerely remain an important and technologically 
relevant branch of research in the near future. 
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